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Abstract. The statistics of the fundamental bulge and disk parameters of galaxies and their relation to the Hubble 
sequence were investigated by an analysis of optical and near-infrared observations of 86 face-on spiral galaxies. The 
availability of near-infrared K passband data made it possible for the first time to trace fundamental parameters 
related to the luminous mass while hardly being hampered by the effects of dust and stellar populations. The observed 
number frequency of galaxies was corrected for selection effects to calculate volume number densities of galaxies with 
respect to their fundamental parameters. The main conclusions of this investigation are: 

1) Freeman's law has to be redefined. There is no single preferred value for the central surface brightnesses of disks 
in galaxies. There is only an upper limit to the central surface brightnesses of disks, while for lower central surface 
brightnesses the number of galaxies per volume element decreases only slowly as function of the central surface 
brightness. 

2) The Hubble sequence type index correlates strongly with the effective surface brightness of the bulge, much better 
than with the bulge-to-disk ratio. 

3) The disk and bulge scalelengths are correlated. 

4) These scalelengths are not correlated with Hubble type. Hubble type is a lengthscale-free parameter and each 
type therefore comes in a range of magnitudes (and presumably a range of total masses). 

5) Low surface brightness spiral galaxies are not a separate class of galaxies. In a number of aspects they are a 
continuation of a trend defined by the high surface brightness galaxies. Low surface brightness galaxies are in general 
of late Hubble type. 

Key words: Galaxies; fundamental parameters - Galaxies: luminosity function - Galaxies: photometry - Galaxies: 
spiral - Galaxies: statistics - Galaxies: structure 



1. Introduction 

The light of a spiral galaxy is dominated by two compo- 
nents, the disk and the bulge. The basic difference between 
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these components lies in their support against gravita- 
tional collapse. The disk is almost completely rotationally 
supported, while the bulge is for some fraction also pres- 
sure supported. At least two parameters are needed to de- 
scribe the light distribution of each of these components: 
a surface brightness term and a spatial scaling factor. 
The fundamental parameters of the disk are usually ex- 
pressed in central surface brightness (/io) and scalelength 
ih), while the bulge parameters are expressed in effective 
surface brightness (/Xe) and effective radius (ro). These 
fundamental parameters were determined for a large sta- 
tistically complete sample of galaxies by de Jong (1995a, 
Paper II). The distributions of the fundamental parame- 
ters are still poorly known and their statistics are investi- 
gated in this paper with some emphasis on three relation- 
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ships: 1) " Freem an's law" , the empirical relation found by 
Freeman ( 1970 ) indicating the constancy of disk central 
surface brightness among galaxies, 2) the number density 
of galaxies as a function of their fundamental parameters 
and 3) the relation between the fundamental parameters 
and Hubble classification. 

1.1. Freeman's law 

One of the most remarkable results presented in the clas- 
sical paper of Freeman ( 1970 ) was the apparent constancy 
of the B passband /iq of spiral galaxies. For a subsample of 
28 (out of 36) galaxies he found (/^o) = 21.65 ± 0.3 B-mag 
arcsec"^. If the central M/L is approximately constant 
among galaxies, this translates directly into a constant 
central surface density of matter associated with the lu- 
minous material. 

Several authors have tried to explain this result. It 
has been argued that ignoring the contribution of the 
bulge t o the light profile could produ ce the effect (Kor- 
mend y |l977| ; PhiUipps & Disney |l983t Davies |199G| ). Free- 
man (1970) did not decompose the luminosity profiles in a 
bulge and disk, but fitted a line to the linear part of the lu- 
minosity profile plotted on a magnitude scale. This linear 
part of the profile could be contaminated by bulge light. 
With their models Kormendy ( 1977 ) and Davies ( 1990 ) 
show that the central surface brightness of low surface 
brightness disks will be overestimated by this procedure 
because of the extra bulge light near the center. The cen- 
tral surface brightnesses of high surface brightness disks 
with a short scalelength are underestimated; because of 
the small disk scalelength the bulge light dominates the 
luminosity profile again in the outer region, but with a 
longer scalelength and a lower surface brightness than the 
disk. Several arguments can be raised against this inter- 
pretation (see also Freeman 1978): 1) even with bulge light 
included the result is still important, 2) many later type 
galaxies hardly have a bulge, but the effect is still present 
(van der Kruit 1987), 3) in samples where proper decom- 
position techniques are used the effect is still found, al- 
though with a larger dispersion (Boroson |l98l| ) , 4) a lim- 
ited range in bulge parameter space was explored in the 
models mentioned above, which might not be representa- 
tive of the bulges in spiral galaxies. 

Dust extinction has also been proposed as an explana- 
tion for the constancy of //q (Jura 198C; Valentijn 1990). If 
galaxies are optically thick in the B passband, one is only 
looking one optical depth into the galaxies and always ob- 
serves the same outer layer. This removes the inclination 
dependence from the Freeman relation, but leaves the un- 
solved problem of why all galaxies should have the same 
surface brightness at optical depth equal to one. 

Freeman established his relation in the B passband 
where the light of galaxies is dominated by a very young 
population of stars, which make up only a few percent of 
the stellar mass. Of all commonly used passbands the light 



of the massive old stellar population is relatively the most 
important in the near-infrared (near-lR) K passband used 
here. The K passband has the additional advantage that 
the extinction by dust is strongly reduced. The K pass- 
band is therefore best suited to trace the fundamental pa- 
rameters of the luminous mass. However, other passbands 
have been used as well in this study to investigate the 
wavelength dependence of the bulge and disk parameters 
due to dust and population effects. 



De Vaucouleurs (1974) was one of the first to suggest 
that the constancy of /ito might result from a selection ef- 
fect. This was later quantified by Disney (1976) and Allen 



& Shu (1979). Catalogs of galaxies have usually been se- 
lected by eye from photographic plates using some kind 
of diameter limit. One might therefore select against very 
compact galaxies with a high central surface brightness, 
because these have small isophotal diameters. Likewise, 
galaxies with a very low surface brightness might have 
been missed due to the lack o f cont rast with the sky back- 
ground. Disney & PhiUipps ( 1983 ; see also Davies 1990 ) 
define a visibility for a galaxy, which enables one to cor- 
rect a sample for these selection effects if one has made a 
careful initial sample selection. 



1.2. Bivariate distributions 

Correcting for selection effects is in fact trying to deter- 
mine from the observed statistics how many galaxies there 
are per unit volume with a certain property. More than 
one property can be used in determining such a distribu- 
tion per volume. One needs at least two parameters to 
characterize the exponential light profile of a disk domi- 
nated galaxy and a bivariate distribution function of both 
disk parameters is a more general statistical description of 
galaxy properties than a one parameter function. The di- 
ameter, the central surface brightness and the luminosity 
distribution functions of galaxies are integrations of this 
bivariate distribution in a certain direction. In this process 
information is lost and the bivariate distribution function 
is therefore more useful in studies of deep galaxy counts 
and provides more constraints on theories of galaxy forma- 
tion and evolution than its one dimensional counterparts. 

Bivariate distribution functions of galaxies have been 



determined only a few times before (Choloniewski |1985 



Phillips & Disn ey |1986| ; van de r Kru it |l987i |1989| ; Saun- 
ders et al. |1990|; Sodre & Lahav|l993D. Even though differ- 



ent fundamental parameters are used, almost all (except 
Saunders et al.) of these distributions describe fundamen- 
tally the same thing in different ways. These studies were 
performed in the B or comparable passbands, which is, 
as mentioned before, not the wavelength most suited to 
study global fundamental properties of galaxies. 
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1.3. Morphological classification 

For classification of spiral galaxies on the Hubble sequence 
three principal discriminators are used: 1) the pitch-angle 
of the spiral arms, 2) the degree of resolution of the arms 
(into Hii regions, dust lanes and resolved stars) and 3) 
the bulge-to-disk (B/D) ratio. In his detailed description 
of the Hubble sequence, Sandage (1961) indicates that the 
B /D ratio is the weakest discriminator unless galaxies are 
seen edge-on. He finds clear mismatches in type between 
classifications using items 1) and 2) and classifications us- 
ing item 3). Another factor hampers the use of B/D ratio 
for classification of early spirals. On the photographs used 
for classification the central region of an early spiral galaxy 
is normally overexposed in order to show clearly the faint 
spiral structure. 

Still, the B /D ratio is often assumed to be the principle 
parameter underlying the Hubble sequence, even though a 
tight correlation between classification and measured B /D 
ratios was never found. The measurements indicate at best 



a trend (e. g. Sim ien & de Vaucouleurs 1986 ; Andredakis 
& Sanders 1994 ) and the discrepancies between B/D ra- 
tio and Hubble type have been attributed to two sources 
of error. First there is the uncertainty in classification. 
Comparisons of Hubble types given by different classifiers 
show an rms uncertainty in type index of order 2 T-units 
(Lahav et al. 1995| ). The second source of error is the un- 
certainty in the bulge/disk decomposition, due to, among 
other things, the mathematical peculiarities of the widely 
used r-^/^ bulge law (de Vaucouleurs 19481 ). 



1.4. Outline 

The main goal of this investigation is to determine the na- 
ture of the Freeman law. In order to address the problems 
concerning the Freeman law, a large sample of face-on spi- 
ral galaxies was carefully selected and surface photometry 
was obtained in the K passband as well as in several other 
passbands. A large number of other global and structural 
parameters of the galaxies were determined in this inves- 
tigation and their nature is also explored in this paper. 

The remainder of this article is organized as follows. 
The data set and the extraction of the observed bulge 
and disk parameters are briefly described in Sect. ^. The 
corrections to the observations in order to calculate num- 
ber distributions are described in Sect. ^ and these distri- 
butions are presented for the B and the K passband in 
Sect. ^. The relations found are discussed within the con- 
text of the three main points of interest (Freeman's law, 
bivariate distributions and Hubble sequence) in Sect. |[ 
The conclusions are summarized in Sect. |^. 

2. The data 

In order to examine the parameters describing the global 
structure of spiral galaxies, 86 systems were observed in 
the B, V, R, I, H and K passbands. A full description of 



the selection, observations and dat a ext raction can be 
found in de Jong & van der Kruit ( 1994 , Paper I). The 
galaxies in this statistically complete sample of undis- 
turbed spirals were selected from the UGC (Nilson |l973| ) 
to have red diameters of at least two arcmin and minor 
over major axis ratios larger than 0.625. The survey was 
limited to 12.5% of the sky globe. Standard reduction 
techniques were used to produce calibrated images. 

In Paper II the extraction of the bulge and disk param- 
eters from the calibrated images is described. An extensive 
error analysis was performed using different fit techniques. 
The best results were obtained with a model galaxy with 
an exponential radial light profile for both bulge and disk, 
that was two-dimensionally (2D) fitted to the full cali- 
brated image. This 2D fit technique made i t also possible 
to fit an additional Freeman bar (Freeman 1966) compo- 
nent, which improved the fit for 23 of the 86 galaxies. 
The error analysis revealed the two dominant sources of 
error in the derived component parameters to be: 1) the 
assumed luminosity profile of the bulge and 2) the uncer- 
tainty in the sky background subtraction. Other uncer- 
tainties, like errors in seeing correction, zero-point errors 
and resolution problems were found to be much smaller in 
most cases. 

Assuming that the exponential profile is a reasonable 
description of the bulge light distribution, the dominant 
source of error in the parameters is caused by the uncer- 
tainty in the sky background level. This uncertainty was 
not taken into account using the 2D fit technique, but 
the ID errors can be used, because the 2D fit results are 
generally comparable to the double exponential ID fit re- 
sults (see Paper II, Fig. 6). The ID errors do include the 
uncertainty of sky background subtraction and are always 
larger than the formal 2D fit errors. The ID errors are only 
shown in the graphs presented here if they are significantly 
larger than the symbol size. 

In this paper the RC3 (de Vaucouleurs et al. |l991 ) 
morphological type index T is used (see also Paper I). Be- 
cause a few galaxies had no RC3 classification, I classified 
them as UGC 1551~(8), UGC 1577-(4), UGC 9024-(8) and 
UGC 10437-(7). The mean error in type index in the RC3 
is stated to be 0.89. This number seems to be very low. 
Lahav et al. (1995) showed that the dispersion between 
the RC3 T-index and the T-values of six expert classifiers 
was on average 2.2 T-units for a sample of 831 galaxies. 
The dispersion between any two classifiers ranged between 
1.3 and 2.3 T-units, with 1.8 on average. It is safe to say 
that the uncertainty in classification in the RC3 is at least 
1.5 T-units. 

The data set comprises 86 galaxies in six passbands. To 
keep a clear view on the obtained results I will concentrate 
on the two most extreme cases, the B and the K passband 
data. The results for the other passbands are available 
in electronic form. The B and K passband results are 
displayed in the graphs with the same dynamic range (but 
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often with different zero-points) and therefore they can be 
compared directly. 

3. Corrections 

The observed bulge and disk parameters determined in 
Paper II have to be corrected for all kinds of systematic ef- 
fects. These corrections are often uncertain but necessary. 
One can only expect that they are at least in a statistical 
sense correct. 

3.1. Galactic foreground extinction 

The measurements of brightness and surface brightness 
were corrected (unless stated otherwise) for Galactic fore- 
ground extinction according to the precepts of Burstein & 
Heiles ( 1984 ) and the actual B passband extinction values 
were adopted from the RC3. The Galactic extinction curve 
of Rieke & Lebofsky ( 1985 ) was used to convert these B 
passband extinction values to other passbands. The sam- 
ple galaxies were selected to have a Galactic latitude larger 
than 25°; the extinction correction is in general small and 
gets smaller for the longer wavelength passbands. The av- 
erage correction is 0.14 i?-mag and the largest correction 
is 0.68 _B-mag, which translates into 0.06 K-vci&g. 

3.2. Inclination corrections 



Since Valentijn ( |199C| ) reopened the debate of optically 
thin versus optically thick spiral galaxies, inclination cor- 
rections for surface brightness have become less trivial. A 
simple equation for correcting surface brightnesses for in- 
clination effects, taking internal extinction into account, 
has the form 



= /i- 2.5Clog(a/5), 



(1) 



where a/b is the major over minor axis ratio of the galaxy 
and C the internal extinction parameter, which takes val- 
ues 0<C<1. Fully transparent galaxies are described by 
C = 1, while the case C = describes the optically thick 
ones. 

It is unlikely that the inclination correction indeed 
takes such a form in the optical passbands, as extinction in 
the optical passbands is for a considerable fraction caused 
by scattering and not just by absorption alone. Light will 
be scattered preferably from edge-on directions to face-on 
directions, which means that extinctions will seem to be 
higher for edge-on than for face-on galaxies. On top of 
that, certain configurations of dust and stars can behave 
optically thin in an inclination test, while they may in fact 
be completely opaque. A clear example of this is a very 
thin layer of optically thick dust between a thicker slab 
of stars. It is not trivial to produce a better description 
as there are too many unknowns and C itself may be a 
function of galactic radius (see e.g. Giovanelli et al. 1994; 
Byun et al. 1994). Therefore Eq. (|l|) is used as a working 



hypothesis. However, for a face-on selected sample such 
corrections are small. The average correction for the sam- 
ple examined here is 0.26 mag arcsec"^ when C = l, with 
a maximum of 0.60 mag arcsec^^ for the galaxy with the 
largest observed h/a of 0.58. 

3.3. Distances 
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Fig. 1. The distance distribution of the sample galaxies. For 
the dashed line the Vqsr velocities from the RC3 were used, the 
full line indicates the distance distribution when the velocities 
are corrected for Virgo-centric infall. 



The distances to the observed galaxies were calculated us- 
ing a Hubble flow with an Hq of 100 km s~^ Mpc~^, cor- 
rected for infall into the Virgo cluster using the 220 model 
of Kraan-Korteweg (1986). This model assumes that the 
Local Group has an infall velocity of 220 km/s towards 
the Virgo cluster and describes the motions of the galax- 
ies around the cluster by a non-linear flow model. The 
Vhei velocities needed for this model were calculated from 
the Vqsk velocities listed in the RC3, which are also tab- 
ulated in Paper I. The nearest galaxy is at 6.2 Mpc, the 
most distant galaxy is at 82.5 Mpc. The peculiar veloci- 
ties of galaxies were assumed to be on average 200 km/s 
in the line of sight, which introduces an uncertainty (ud) 
of 2 Mpc in the distance estimates. The distribution of 
distances is displayed in Fig. |l|, which shows some excess 
of galaxies at ~45 Mpc because of an extension of the 
Pisces-Perseus supercluster. The relationships presented 
in this study are very little affected when other infall and 
flow models are used to calculate distances. 



3.4. Selection correction 

The physically relevant quantities are not the observed 
numbers of galaxies with a certain property, but the fre- 
quency of galaxies with a certain property in a volume. 
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Therefore, the fact that a galaxy is included in the sam- 
ple has to be linked to the statistical probability of find- 
ing such a galaxy in a certain volume. The galaxies in 
the sample were selected to have UGC red major axis 
diameter {D]^^) of at least 2 arcmin. This creates a se- 
lection bias against galaxies with low surface brightness 
and/or small scalelengths, as they appear smaller on pho- 
tographic plates. The distances (d) to the observed galax- 
ies and their angular diameters (I?maj) are known and 
the maximum distance at which galaxy can be placed, 
while still obeying the selection criteria, can be calculated 
(rfmax = rf£'maj/-D|™j). A galaxy can only enter the sam- 
ple if it lies in a spherical volume which has this max- 
imum observable distance as radius. Turning this argu- 
ment around, one can expect on statistical grounds that a 
selected galaxy samples a spherical volume with a radius 
equal to its maximum observable distance (a more for- 



mal discussion can be found in Felten 1976). The volume 
sampled by a galaxy in a diameter limited sample is thus 



An 

n 



47r 



(2) 



Following the previous line of reasoning, an estimate for 
the average number of galaxies in a unit volume obeying 
a certain specification (S) for a complete sample of N 
galaxies is 



N 

E 



(3) 



where i is summed over all N galaxies in the sample and 
— 1 if the specification is true for galaxy i and — 
if false. The error in ^{S), assuming Poison statistics in 
a homogeneous universe and considering the uncertainties 
in the distances, can be calculated by 



N 

E 



N 



(4) 



There is always a chance that a member of a peculiar 
class of galaxy happens to be nearby and gets a lot weight 
in Eq.(|3|) and this volume correction can therefore only 
be applied to large samples. One must ensure that a large 
enough volume of space is sampled so that galaxies are 
randomly distributed in space. Figure |l| shows that the 
sample mainly traces the local density enhancement, as 
large scale structures in the universe have scales of order 
50 Mpc. Equation H should therefore be used with care, 
because the number of galaxies with small intrinsic diam- 
eters will be overestimated relative to the larger ones due 
to the local density enhancement. The average number of 
galaxies per Mpc'^ calculated with Eq. (||) might be more 
representative of the local environment than of the mean 
cosmological values. Still it is a useful equation to observe 
general trends in bivariate distributions and to compare 
results obtained from different passbands. 



Other methods to correct distributions for selection ef- 
fects have been advocated, because they take spatial den- 
sity fluctuations into account (for an overview see Efs- 
tathiou et al. 198^). These methods assume that the in- 



trinsic distribution function is independent of position (x) 
in space, so that we can write $(S') = (j>{S)p{'x.), thereby 
losing the absolute calibration of the number density. 
These methods all assume a clear relation between the dis- 
tribution parameter(s) and the limiting selection param- 
eter(s). This is not the case for the current investigation. 
A diameter limit is not trivially linked to the central sur- 
face brightness distribution, certainly not when a different 
passband is used for the selection and the distribution. 

The correction of Eq. (^) is only valid if a particular 
galaxy would have been measured at the same intrinsic (as 
opposed to angular) diameter, had it been at a different 
distance. In Paper I it was shown that this is probably the 
case for the UGC galaxies with type index T<6. For later 
types the situation is less clear, there is a too short a range 
in diameters to check and it must be assumed that for 
late-type systems the same type of galaxy is measured at 
the same intrinsic diameter at different distances. Under 
this assumption it is not important that the UGC eye 
estimated diameters of late-type galaxies correspond to 
lower average surface brightness than that of early types 
(see Paper I, Fig. 11). This effect just means that there are 
more late-type galaxies in the sample than expected based 
on their isophotal diameter, but their average distance 
will be larger so that the number of galaxies per sampled 
volume stays the same. 

The volume correction of Eq. (||) can be used to cal- 
culate number density distributions for all passbands, as 
long as the red UGC diameters are used to calculate the 
Vraa,^. The distribution of any galaxy parameter 5' can be 
determined in any passband; the use of the red UGC di- 
ameters in Eq. (|^) ensures the correction for the intrinsic 
selection effects of the whole sample. 

Next to the diameter limit, there are two more selec- 
tion criteria defining the sample. The selection was limited 
to 12.5% of the sky and only galaxies with 5/a> 0.625 were 
used, which is only 37.5% of all possible random orienta- 
tions. Equation ^ was corrected for these selection criteria. 
A correction was also applied for the fraction of galaxies 
for which no (photometric) data was available in a certain 
passband. All these corrections were made under the as- 
sumption that the incompleteness had no correlation with 
the investigated parameters. 

Equation (|^) can only be applied when the sample is 
complete. The statistical completeness of the sample can 
be tested with the V/Vma^-test (Paper I). The V/Vmax of 
a galaxy is the spherical volume associated with the dis- 
tance of a galaxy divided by T4iax as defined in Eq.(||), 
thus for a galaxy in this diameter limited sample V/Vmax = 
(-^maj/^™'^j)'^' objects distributed randomly in space 
the average value of T^/Knax should be 0.5 ± 1/V12 x N, 
where N is the number of objects in the test. For the cur- 
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rent sample (F/Vmax) = 0.57 ± 0.03 and therefore there 
are slightly too many galaxies with a small angular di- 
ameter in the sample. The original sample of 368 galax- 
ies from which the current subsample was selected had a 
(^/Kiax) = 0.496 ± 0.015 (Paper 1). Subsequent selection 
depended only on the position on the sky and therefore 
the excess of small diameter galaxies is probably caused 
by the density enhancement of the Pisces-Perseus super- 
cluster, which gives some extra galaxies at the diameter 
selection limit. This might give some extra high surface 
brightness and/or large scalelength galaxies in the sample 
above the cosmological mean, because galaxies have to be 
intrinsically large to be included in the sample being at 
the distance of the Pisces-Perseus supercluster. 



In a recent paper Davies et al. (1994) argued that the 
sample used by van der Kruit (1987) was incomplete in a 
magnitude y/Knax-test. They argued that a hidden mag- 
nitude limit had influenced the selection, so that an extra 
selection correction should be applied. I will follow up on 
this argument as the sample used here has been selected 
Qsing similar criteria as van der Kruit used for his sample. 

There is nothing hidden about a magnitude selection 
effect for a diameter limited sample. On the contrary, it 
is expected. For galaxies with a certain absolute magni- 
tude M there exists a range of possible {iiQ,h) combi- 
nations, satisfying Af cx /xq — 5 log(/i), but only a limited 
range of them will satisfy the diameter selection criterion 
-DmajCx(Miim^A'o)'i>^maj- Thus for galaxics of the same 
apparent magnitude we will miss some of the small scale- 
length, bright /io and some of the large scalelength, faint 
fiQ galaxies, while still having selected a complete sample 
in diameter. The complication arises because apparent di- 
ameters and magnitudes are not independent parameters 
for galaxies and their V^/Vi„ax-tests cannot be applied in- 
dependently. Similar to Eq. a Vmax corresponding to 
a magnitude limit can be constructed. The smallest of the 
magnitude V^ax and the diameter V^ax values should be 
used for each galaxy in a combined F/T^nax-test. These 
tests cannot be performed separately. The fundamental 
premise of a diameter limited sample is that the diame- 
ter Vmax is always smaller than the magnitude Vmax, and 
therefore a "hidden" magnitude limit does not have to be 
taken into account and Eq. (H) is sufficient. 



Davies et al. ( 1994 ) also indicate that van der Kruit's 
sample becomes incomplete for low surface brightnesses at 
Ho > 22.3 as (F/Kaax) = 0.35 ± 0.08. I would like to note 
that this might just be a statistical fluctuation of low num- 
ber statistics, as (F/Kiax) = 0.41 ± 0.10 for /io > 22.5 and 
= 0.46 ± 0.13 for /^o > 22.7, and thus for even lower 
surface brightnesses the sample is in the statistically complete 
range of {V/Vmax> = 0.5. 
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Fig. 2. The apparent scalelength versus observed central sur- 
face brightness of the disks. The dashed line indicates the se- 
lection limit of Dmaj > 2' at the 24.7-R-mag arcsec"^ isophote 
for galaxies with perfect exponential disks. The indicated error 
estimates are the errors from the ID profile fitting with expo- 
nential bulge taking the uncertainty in the sky background 
level into account. The errors are only plotted if they are sig- 
nificantly larger than the symbol size. 



4. The distribution of disk, bulge and bar param- 
eters 

In this section 1 investigate the distributions of the struc- 
tural parameters of the different galaxy components as a 
function of morphological type and of each other. First the 
structural parameters of the disk and bulge are examined 
independently. In the final subsection, the relationships 
between disk and bulge parameters are investigated. The 
distributions of bulge and disk parameters are corrected 
for selection effects to yield volume number densities. 

4-.1- The disk parameters 

Figure ^ indicates some aspects of the completeness and 
selection effects of the sample. It shows the distribution 
of the observed central surface brightnesses versus scale- 
length as obtained from the 2D fits of Paper 11. The R 
passband values are plotted, because these values are most 
closely related to the (red UGC diameter) selection cri- 
teria. The dotted line indicates the selection limit for a 
diameter cutoff at 2 arcmin at a surface brightness of 
24.7 i?-mag arcsec"^ for a perfect exponential disk. The 
24.7 i?-mag arcsec"^ is the average surface brightness at 
which the UGC red diameters were determined (see Pa- 
per I). As mentioned in Paper I, not all UGC galaxies had 
their diameters estimated at the same isophote level. This 
explains why there are some galaxies to the left of the 
selection line in Fig. ||. If all galaxies have the same scale- 
length, the number of galaxies expected in the sample will 
decrease as h'^^ and therefore it is not surprising that there 
are hardly any galaxies in the sample below 22 i?-mag 
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(^o) (-B-mag arcsec 


') 




(mo> 


(if- mag arcsec 


') 


RC3 type 


nr. 


C = 


C = 0.5 


C=l 


nr. 


C = 


C = 0.5 


C = 1 


< T < 6 


61 


21.32 ± 0.78 


21.45 ± 0.76 


21.58 ± 0.74 


60 


17.48 ± 0.71 


17.61 ± 0.69 


17.75 ± 0.67 


6 < T < 8 


12 


22.01 ± 0.75 


22.16 ± 0.73 


22.30 ± 0.72 


10 


18.34 ± 0.90 


18.50 ± 0.90 


18.65 ± 0.90 


8 < T < 10 


8 


22.97 ± 0.60 


23.12 ± 0.57 


23.26 ± 0.55 


7 


20.05 ± 1.05 


20.21 ± 1.03 


20.37 ± 1.01 


all 


81 


21.59 ± 0.92 


21.72 ± 0.90 


21.86 ± 0.89 


77 


17.82 ± 1.08 


17.96 ± 1.08 


18.10 ± 1.08 



Table 1. The average Galactic extinction corrected central surface brightnesses for different inclination corrections (Eq. l)^) 
and type index bins. C = Q corresponds to an optically thick disk, C — 0.5 to a semi transparent disk and C = 1 to a fully 
transparent disk. The values are in mag arcsec"^ with their standard deviations. 
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Fig. 3. The Galactic extinction corrected central surface brightness of the disks as function of morphological RC3 
crosses show the values averaged over the bins indicated by the horizontal bars. The vertical bars indicate the standard 
of the mean values. 
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arcsec ^ . Obviously no galaxies can enter the sample with 
^0 fainter than ^^24. 7 i?-mag arcsec"^. 

Let us now look at the central surface brightness as 
function of morphological type (Fig. Apparently the 
galaxies from type T = 1 to 6 have on average the same fiQ , 
but with a large scatter. The later types have on average 
a significantly lower central surface brightnesses and they 



might be classified as late types just because they are low 
surface brightness (LSB) systems. It can be readily seen 
that this difference between early and late-type galaxies 
increases when going from the B to the K passband. This 
indicates that disks of the later type spirals are bluer than 
the disks of the early ones, but the discussion on the colors 
of these galaxies is postponed to Paper IV in this series 
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(de Jong 19951: ) . The average /io values were calculated for 
three morphological type bins indicated by the horizontal 
bars in Fig. ^ as well as for the total sample. The values 
with their standard deviations are tabulated in Table ^ 

The average /io values were also calculated with an in- 
clination correction according to Eq. (|l|) with values for 
C = 0.5 and C — 1 (semi transparent and completely 
transparent behavior). The results can also be found in 
Table |l|. The standard deviations on the average /io val- 
ues are slightly smaller for C = 1 , and even though it is a 
small effect, it is persistent for all subgroups and all pass- 
bands. The main result is of course a shift in the mean 
central surface brightness of the disks. For all remaining 
plots an inclination correction with C= 1 will be used. 

The distribution of the other disk parameter, the scale- 
length (h), as function of type is shown in Fig. ^. There is 
no trend of h with type and there is a large range in scale- 
lengths. There might be a lack of late-type galaxies with 
small scalelengths, but this can probably be attributed 
to a selection effect: the selection criteria are heavily bi- 
ased against LSB galaxies with small scalelengths. The 
scalelengths are smaller in the B passband than in the K 
passband (discussion in Paper IV) . 

The information of Figs ^ and ^ are combined in Fig. |[ 
This figure shows that there is an upper limit in the (/io,/i)- 
plane, as there are no galaxies with large scalelengths and 
high central surface brightnesses. This cannot be caused 
by selection effects, large bright galaxies just cannot be 
missed in a diameter selected sample. This upper limit 



has been noted before by Grosb0l (1985). The upper limit 
partly follows the line of constant total disk luminosity, as 
indicated by the dash ed line in Fig. ||. Note that the TuUy- 
Fisher relation ( 1977 , hereafter TF-relation) implies that 
this is also a line of constant maximum rotation speed 
of the disk. There is also an upper limit to the central 
surface brightness at about 20 i?-mag arcsec^^ (16 K- 
mag arcsec~^). Again galaxies brighter than these limits 
are hard to miss because of selection effects. 

Late-type galaxies have lower central surface bright- 
nesses in Fig. 1^, but the early and intermediate types show 
no segregation. The scalelengths also gives no segregation 
according to type. Very few late-type galaxies with very 
short scalelengths were selected, but as shown before, late- 
type galaxies have lower surface brightnesses and the se- 
lection biases against galaxies with low surface brightness 
and short scalelengths are large. These biases are indicated 
by the dotted lines in Fig. ||. To the right of these lines the 
sample should be complete to the indicated distance. The 
lines are calculated under the assumption that all galaxies 
have perfect exponential disks with the same color at the 
selection radius {B-R=1.3, R-K = 2.5) and that the se- 
lection limit is at 2' diameter at the 24.7 i?-mag arcsec"^ 
isophote (as in Fig. ||). Although these assumptions are 
not valid for an individual galaxy, the dotted lines help to 
estimate the selection effects; the galaxies near the 50 Mpc 
line had about 125 times more chance of being included 



in the sample than the galaxies near the 10 Mpc line! The 
fact that the number density of objects does not decrease 
by 125 from one line to the other already indicates that 
there are many more "small" galaxies per volume element 
than "large" galaxies. 

The distribution of the absolute magnitude of the disk 
(Mdisk) against type (Fig.^) can also be deduced from 
Figs ^ and ^ (Afdisk c>c /io — 2.5 log(27r/i^), no inclination 
dependent extinction correction was applied). As scale- 
lengths show little correlation with type, the distribution 
of disk magnitudes reflects the distribution of the central 
surface brightness against type. There was no apparent 
segregation according to bar classification in Figs ^, ^ || 
and |g. 

So far, only the observed distributions were presented, 
but the distributions per volume are of more importance. 
Therefore the volume correction as described in Sect. ^ 
was applied. The correction transforms Fig. ^ into the bi- 
variate distribution in the {^Q,h)-planc presented in Fig. ^ 
This is a representation of the true number distribution 
of spiral galaxies per volume element of one Mpc'^ with 
respect to both disk parameters. The magnitude and /io 
upper limits noticed in Fig. ^are also present here. We are 
dealing with low number statistics now, which is reflected 
in the erratic behavior of the distribution. The uncertainty 
increases in the direction of small scalelength and low sur- 
face brightness. These galaxies have so small isophotal di- 
ameters that they really have to be nearby to be included 
in the sample and such a small volume is sampled that 
statistics are working against us. For example if the true 
volume densities in the (17i^-mag arcsec"^, 1 kpc) and 
(21 if -mag arcsec^^, lkpc)-bins are equal, the chance of 
observing a galaxy in the last bin would be 0.5. If there 
had been such a galaxy in the sample, a lot of weight 
would have been given to it. In short, the distributions 
are not well sampled in the low surface brightness, small 
scalelength region. No galaxies were selected in this region, 
but the traced volume is also very small. The dominant 
type of spiral galaxy has a scalelength of about 1 kpc and 
a central surface brightness of 21 i?-mag arcsec"^ (17 K- 
mag arcsec"^). 

By summing all bins in one direction, the bivariate 
distributions of Fig. ^ can be used to calculate the dis- 
tributions of fiQ and h separately. This figure indicates 
therefore where one can expect problems in the determi- 
nations of the /io and h distributions due to the under- 
sampling in the low surface brightness, small scalelength 
region. The /io distributions will get incomplete for cen- 
tral surface brightnesses fainter than 21.5 B-mag arcsec"^ 
(19 K-ma.g arcsec"^) and the h distributions should not 
be trusted for scalelengths smaller than 1 kpc. The under- 
sampling in the fio distribution is considerably reduced 
when only the galaxies with scalelength larger than 1 kpc 
are used. The undersampling problem of this sample could 
also be circumvented by imposing absolute magnitude or 
intrinsic diameter limits. To speak of the /io distribution is 
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Fig. 6. The distribution of absolute disk magnitudes as function of type index. 
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incorrect and one should indicate to what type of galaxies 
the sample is restricted. 

The distributions of central surface brightnesses of 
galaxies with scalelength larger than 1 kpc are displayed in 
Fig. ^. The distributions are remarkably flat for the total 
sample. The number density density decreases by about a 
factor of 4 from /io* — 21 to 24 i?-mag arcsec"^ and by a 
factor ~10 from 17.5 to 22 K-mag arcsec"^. The distribu- 
tions are narrower when only types earlier than T=6 are 
used. Disks of late-type galaxies are bluer, which makes 
the overall distribution narrower in B than in K . The dis- 
tributions are not limited by selection effects at the bright 
end, even if one assumes there is an upper limit to the 
total luminosity of a galaxy (see Fig. ||). The number den- 
sity of galaxies decreases sharply with /Iq brighter than 20 
B-mag arcsec"^ (~16 K-mag arcsec"^). At the faint end 
a limited volume is sampled, and Fig. indicates that the 



sample is biased against galaxies with a /Iq fainter than 23 
B-mag arcsec"^ even for galaxies with scalelength larger 
than 1 kpc. Obviously galaxies with central surface bright- 
ness fainter than 26 i?-mag arcsec"^ could never enter the 
sample. The distributions of jiQ in Fig. ^ could be slightly 
higher at the faint end and should probably be extended 
to much lower surface brightnesses. 

The volume corrected distributions of the logarithm 
of the scalelengths (Fig. ^ show first a small increase of 
galaxies to scalelengths of about 1 kpc. This is probably 
caused by the undersampling effect at low surface bright- 
nesses and small scalelengths. For scalelengths larger than 
1 kpc we notice a steady decline of about a factor 100 in 
one dex. There is no segregation with morphological type. 

The most important results obtained in this subsection 
are as follows. There is a large range in disk central sur- 
face brightnesses among galaxies, mainly due to the lower 
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Fig. 10. The Galactic extinction corrected effective surface brightness of the bulge as function of morphological RC3 type. 
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Fig. 11. The effective radius of the bulge as function of morphological type. 
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surface brightnesses of late-type galaxies. The range de- 
creases slightly when a transparent inchnation correction 
is used. Selection effects are very significant in determin- 
ing number density distributions and after correcting for 
these effects there is no single preferred value for the cen- 
tral surface brightnesses of disks. There is an upper limit 
to /io, but the number density distribution decreases only 
slowly at the faint end. 

4-2. The bulge parameters 

The same diagrams used to describe the disk parameters 
are now used to present the bulge parameters. The dis- 
tributions of effective surface brightness are presented in 
Fig. |l^. The effective surface brightness shows a tight cor- 
relation with type index, especially considering the uncer- 
tainty of at least 1.5 T- units (1 sigma) in type index (La- 
hav 1995 ). Almost all of the scatter can be explained by 
this uncertainty. This relation also holds for the /ic param- 
eters obtained with the other fitting methods presented in 



Paper II, although with a slightly larger scatter. There 
is no apparent correlation of effective radius with galaxy 
type (Fig. |ri|). The relations in Figs ^ and ^ are tighter 
in K than in B. There are several explanations for this 
effect. Bulges are relatively brighter with respect to the 
disks in K compared to B, which will make the fit routine 
work better. Furthermore, circumnuclear star formation 
and dust lanes will affect the B passband more than the 
K passband and make the quality of the decomposition 
worse. There are some galaxies in the sample with clear 
circumnuclear star formation and with dust lanes right 
down to the center. Finally, the re is t he effect of the Free- 
man Type II profiles (Freeman 197C) which is reduced in 
K, thus making fitting easier (see Paper II). The distri- 
bution of points in the (/ic,re)-plane (Fig. |l^) shows no 
correlation. The absence of a correlation between rc and 
morphological type makes the trend in the distribution of 
the absolute bulge magnitude (Mbuige oc /Zc — 2.5 log(r^)) 
versus type (Fig. O) dominated by the Hc- 
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The bivariate distribution of and (Fig. ^4|) shows 
no trends. The dominant type of galaxy in our local uni- 
verse has a bulge with effective radius in the range of 0.1- 
0.3 kpc and effective surface brightness of order 21 B-mag 
arcsec"^ (~16 i^-mag arcsec"^). The relation between the 
bulge parameters and the diameter selection criterion is 
not very obvious and therefore all galaxies are used in the 
calculations of the separate bulge parameter distributions. 
The volume corrected distributions of the bulge fi^ and 
(Figs |l^ and |l^) show the same behavior as the disk pa- 
rameters, i.e. constant distribution of the effective surface 
brightness and a steady decline of a factor of 50-100 per 
dex of the effective radius in the K passband. 

4-3. The bulge/disk relation 

The chronology of the bulge and disk formation is a major 
issue and the relationships between bulge and disk param- 



eters might give some insight in this matter. A strong cor- 
relation between bulge and disk parameters is expected 
if the bulge formed from the disk by secular evolution. 
A correlation might be expected in the hierarchical infall 
and small merger models producing bulges, because both 
bulge and disk originate from the same smaller compo- 
nents. In models where the bulge forms first and the disk 
forms later, there is no obvious reason for a bulge-disk 
correlation. 

Comparing the fj,Q with (Fig. |l^) we see no corre- 
lation, except that the late-type spirals clearly stand out. 
This is most obvious in the K passband. The correlation 
between h and (Fig. |l^) is only tight in the K pass- 
band and not in the B passband. Actually the correlation 
is becoming steadily tighter from the B to the K passband 
with correlation coefficients increasing from 0.6 in B and 
to ~ 0.75 in R and / and ~ 0.8 in H and K passbands. 
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Fig. 17. The central surface brightness of the disk versus the effective surface brightness of the bulge. Different symbols are 
used to denote the indicated morphological type ranges. 
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Fig. 20. The bulge to disk ratio as function of morphological type, using the results from the 2D fit of Paper II. 
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Fig. 21. The bulge to disk ratio as function of morphological type using an r^^^ law bulge in the bulge/disk decomposition. 



The equation for the least squares fitted line is in the K 
passband 

log(rf ) = 0.95 log(/i^) - 0.86 (5) 

with a standard deviation of 0.17. The scalelength dif- 
ference between bulge and disk thus is of order 10. The 
relation also holds for all the ID fit techniques presented 
in Paper II (but less strongly) except for the case of an 
r^/'* law bulge. In the case of ID r^/^ law bulge there is 
at best a weak trend (correlation coefficient 0.15). This 
relation could partly be produced by the fitting routine if 
the errors in both parameters are intrinsically correlated. 
The facts that this relation holds in the K but not in 
the B passband and for both of the totally different ID 
and 2D fitting techniques indicate that the correlation is 
intrinsic and not an artifact of the fit routines. The 
distribution around the solutions found by the fit routine 
also showed no correlation with the relationship between 
disk scalelength and bulge effective radius. 



The absolute magnitudes of bulge and disk correlate 
well (Fig. |l^. This is probably an example that large 
galaxies have more of everything, more bulge, more disk. 
The combined effect of good correlation between r^. and 
h and weak correlation between fj,o and He ensures the 
correlation between Md and Mb. Looking at the bulge- 
to-disk ratio as function of morphological type (Fig. pO| ), 
one sees that there is a correlation, but this correlation 
is less strong than for instance the one of /io with type 
(Fig. The B/D ratio is on average higher in K than 
in B; bulges are redder than disks. This partly explains 
why the correlations which involve bulge parameters are 
tighter in K. The bulge/disk decomposition is more easily 
performed when the bulge is relatively brighter. 

Fig. ^ shows also that the selected galaxies are indeed 
disk dominated systems. The B/D ratios plotted here are 
much smaller than the ratios normally found in the lit- 
erature. This is mainly due to the use of an exponential 
bulge. In Fig. |l] the results with an r^/^ law bulge (Pa- 
per II) are shown. The B/D ratios are higher, but the 
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scatter has increased and there still is no tight correlation 
with morphological type. 

5. Discussion 

In this section I will place the previously described rela- 
tions in the context of the three topics of main interest: 
1) Freeman's law, 2) bivariate distributions and 3) the 
relation between Hubble classification and the structural 
parameters. 

5.1. Freeman's law 



Since Freeman 1970 found disk central surface bright- 



nesses to be constant among spiral galaxies, a number 
of explanations have been brought forward. In the intro- 
duction three explanations were mentioned: 1) optically 
thick dust, 2) erroneous measuring of the disk parameters 
from the light profiles and 3) selection effects. For each of 
these possibilities I check if they are of importance for the 
current sample and whether they can explain Freeman's 
result. 

5.1.1. Optically thick dust 

It has been suggested that opticall y thic k dust co uld be 
the cause of Freem an's law (Jura 1980 ; Valentijn 1990| ; 
Peletier et al. 1994). This is only a partial explanation, 
because it removes the inclination dependence from the 
law. To produce Freeman's law in this way, galaxies must 
have the same surface brightness at r = 1 (where the typ- 
ical surface brightness is produced in an optically thick 
system), which means the problem is only shifted from 
one part of the galaxy to another. One now has to explain 
why all galaxies have the same surface brightness at t — \. 
Taking a constant dust-to-stellar light ratio will not pro- 
duce a constant surface brightness. This is only the case 
if all dust is in front of the star light, but dust and stars 
are mixed in a galaxy and a fraction of stars to the near 
side is less obscured. The amount of extinction in a galaxy 
is not linearly dependent on the amount of dust present 



(de Jong 1995c). Coupling the amount of dust and stars in 
galaxies can only reduce the scatter in /xq, but can never 
produce a constant /xq. 

I have shown that applying the inclination correction 
of Eq. (jl|) reduces the scatter in the /io of the disk go- 
ing from C = to C = 1 (Table |^), indicating transpar- 
ent behavior. The effect is small and the scatter is still 
dominated by the intrinsic differences in the brightnesses 
of the disks. One should realize that the disk parameters 
are largely determined by the outer regions of the galaxy. 
They probably do not refiect the optical thickness of the 
central regions. I note again that galaxies can behave op- 
tically thin in an inclination test, while in fact being op- 
tically thick. 

The K passband data should hardly be affected by 
dust extinction. Looking at Table |l| one can see that the 



standard deviation of /io is for the early-type galaxies 
smaller in the K passband than in the B passband, con- 
trary to what is expected for optically thick dust. The 
increase in standard deviations for the later types can be 
explained by stellar population differences (Paper IV) . I 
conclude that dust extinction is not a major effect in Free- 
man's law, certainly not in the K passband data used here. 

5.1.2. Erroneous profile fitting 



Kormendy (1977), Phillipps & Disney (1983) and Davies 
(1990) have argued that Freeman's law results from fitting 
exponential disks to light profiles without taking the bulge 
contribution to the profiles into account. To prove their 
point, they created model profiles with r^^^ law bulges 
and exponential disks with a range of properties, to which 
exponential disks were fitted in a specified range. These 
models were able to reproduce Freeman's "universal" cen- 
tral surface brightness value of 21.65 i?-mag arcsec^^ with 
just a small scatter. 

The parameter space explored in the models is not very 
physical according to current insights. Kormendy (1977) 
used B/D ratios of 1-120 for the low surface bright- 
ness systems and Phillipps & Disney ( |1983[ ) assumed that 
bulges were so extended that they dominated the light 
profiles again at the 24.5 _B-mag arcsec"^. Davies (1990) 
used a range of properties for the bulge parameters which 
are typical in samples of galaxies, which suffer from severe 
selection effects. He used a constant central surface bright- 
ness of the bulge to show that the central surface bright- 
nesses of disks need not be constant. A change in Vc/h 
ratio was used to produce a range in bulge-to-total light 
ratios {BT). I have shown that /ic is not constant (Fig. 10) 
and that the Vc/h ratio is nearly constant (Fig. ^8|). Even 
though these results where obtained with a exponential 
rather then a r^^'^ law bulge, a constant central bright- 
ness for bulges is excluded and a relationship between re 
and h might be expected. 



I use the parameterization of Davies (199C) to show 
that erroneous bulge/disk decomposition is not a major 
factor in the Freeman law. Figure ^ was produced in the 
same way as Davies' Fig. 6 by fitting exponential profiles 
in the range of 22 to 25.5 mag arcsec^^ to model pro- 
files, which were a combination of an r^/* law bulge and 
an exponential disk. In this figure the intrinsic /iq of the 
model profile is compared with the central surface bright- 
ness of the fitted disk [ijLx)- The Tc/h ratio was taken fixed 
to 0.4 and /ic was adjusted to produce BT ratios in the 
range from 0.05 to 0.75 in steps of 0.1 (contrary to Davies, 
who used a fixed /^c and varied the r^/h ratio). The ten- 
dency to shift intrinsic bright disks to the observed value 
of 21.65 mag arcsec"^ has disappeared. The low surface 
brightness disks have too bright jix values for their /xq, but 
these were clearly fitted in the curved part of the profile 
at the brighter end. From Fig. ^ one can conclude that 
it is unlikely that the high surface brightness disks were 
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underestimated (even with using the "marking the disk" 
fit of Paper II). The situation for LSB systems is not as 
bad as it seems, because the B/D ratios are low for LSB 
galaxies and the curvature of the profiles can be readily 
seen. Choosing values for Vc/h in the range from 0.1 to 1 
hardly changes these results. 



X 

=1 




Fig. 22. The variation of the measured extrapolated central 
surface brightness {^ix) with bulge-to-total light ratio (BT) 
ranging from 0.05 to 0.75, at a given model central surface 
brightness ixq (see text). 



The central surface brightnesses of the galaxies were 
determined in Paper II with a full 2D decomposition tech- 
nique, but also with the "marking the disk" technique and 
a few ID decomposition techniques. If Freeman's result 
was caused by his use of the "marking the disk" method, 
the central surface brightnesses obtained with this method 
should show large and systematic differences with the re- 
sults of the other methods. In Paper II it was shown that 
making a proper decomposition of the profiles in a r^^^ 
law bulge and an exponential disk adds a scatter of at 
most 0.4 mag arcsec"^ to /io with respect to the "mark- 
ing the disk" results. Assuming that about the same value 
would hold for Freeman's sample, this would still result in 
a rather small range in /io for his sample. In Paper II it is 
furthermore argued that the 2D decomposition technique, 
using exponential light profiles for both disk and bulge, 
is more accurate and reduces the rms difference between 
the "marking the disk" method and the 2D fit to 0.3 mag 
arcsec"^. Both the model decompositions using Davies' 
method and the comparisons between different decompo- 
sition methods of real galaxies indicate that it is unlikely 
that Freeman's results were caused by improper decom- 
positions. 

5.1.3. Selection effects 

After taking selection effects into account. Fig. || shows 
that there is no such thing as a simple Freeman's law for 



galaxies with scalelengths larger than 1 kpc. There seems 
to be a clear upper limit to the central surface bright- 
ness, which cannot be explained by selection effects. Even 
taking the apparent upper limit in absolute luminosity in 
Fig. H into account, there still should have been galax- 
ies with /io brighter than 20 B-mag arcsec"^ (16 X-mag 
arcsec"^) in the sample according to this figure. Figure || 
shows further that there is no strong decrease in the num- 
ber of galaxies with lower surface brightness. The distribu- 
tion becomes narrower if we exclude late- type spirals, but 
this exclusion can hardly be justified. Late-type spirals are 
in many respects no separate class of galaxies, but just a 
continuation of the trends set by the earlier type spirals. 
A clear example of such a trend is seen in Fig. 
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In Paper II it was shown that there is at most ~0.30 
mag arcsec"^ rms uncertainty in the central surface 
brightnesses. The uncertainties also showed no correlation 
with the surface brightness itself (Paper II, Fig. 4) and the 
results presented here can therefore not be the results of 
measurement errors. The central surface brightness distri- 
bution of Fig. H changes in some details if one of the other 
fit techniques of Paper II is used, but the general trend re- 
mains unchanged. The same holds true when the distances 
of the galaxies are calculated with other flow models. 

Sample selections are influenced by both /to and h 
and therefore the most important distribution for disk 
dominated galaxies is the bivariate distribution in the 
(/<o,/i)-plane (Fig. 0). These two parameters describe a 
large fraction of the light of disk dominated galaxies and 
to derive this distribution one needs distances. To derive 
the distribution of /io of sample of galaxies without know- 
ing the distance to the galaxies, one must assume that the 
distribution of /io is not correlated to for instance h ov M 
(Davies et al. [l994[ McGaugh et al. |1995| ). The total //o 
distribution can only be calculated in this statistical way 
if the /^o distribution at each h or at each M has the same 
shape. Figure]^ shows that this is probably not the case 
for h and Figure |2^ shows the same for M. The statistical 
methods can at best only be used to get an impression of 
the /io distribution. 

5.2. Bivariate distributions 

The reason for our limited knowledge of low surface bright- 
ness galaxies is clearly indicated by the selection lines in 
Fig. Ij. The use of catalogs like the UGC prevents galaxies 
with central surface brightness fainter than ^25 _B-mag 
arcsec"^ from being included in a sample, independent 
whether the sample is diameter or magnitude selected. 
Only the use of deeper photographic plates enabling a 
selection at fainter isophotes (Schombert et al. 1992| ) or 
deep CCD surveys will result in samples with a larger 
number of low surface brightness galaxies. A galaxy like 
Malin I (Impey & Bothun |l98|) , with /io ~ 26.5 B-mag 



arcsec and ft. ~ 55 kpc, is not found in conventional cat- 
alogs, even though it has an integrated magnitude com- 
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parable to M 101 and a ~10 times as a large scalelength 
as M 101. Figure indicates that galaxies like Malin I are 
probably not very numerous in the local universe, but this 
cannot be said of galaxies with /ip > 24 _B-mag arcsec"^ 
and /iw 1 kpc. There is a clear need for deeper local sur- 
veys, especially in the near-IR. 

Few bivariate distributions have appeared in the liter- 
ature which can be used in comparisons with the current 



results. Van der Kruit (1987) calculated a bivariate distri- 
bution of spiral galaxies in the (fiQ-h) -plane. His distribu- 
tion shows similar features as the distribution presented 
here. The distribution has an upper limit in central sur- 
face brightness at about 21 J-mag arcsec"^ (photographic 
J passband, which is similar to the Johnson B passband) 
and an exponentially declining density distribution with 
scalelength. 

Another comparison can be made with the bivariate 
distribution of van der Kruit ( |1989| ). This distribution 
was constructed in exactly the same way as Fig. |^, with 
only a modification to Eq. (||) to include the effects of an 
additional magnitude selection limit. The central surface 
brightnesses were convolved with a Gaussian of 0.3 mag 
arcsec"^ to incorporate the effects of calibration uncer- 
tainty. This explains the smoothness of the distribution. 
The galaxies of type later than T = 5 were excluded and 
a Hubble constant of 75 km s^^ Mpc~^ was used. The 
same upper limits in absolute magnitude and central sur- 
face brightness can be seen as in Fig. |^. There is only 
one galaxy with /zq brighter than 18 i?-mag arcsec"^ at 
0.7(100/75) kpc. The general trend in the bivariate distri- 
bution of van der Kruit agrees quite well with Fig. ^ with 
equal probabilities along lines of equal absolute luminos- 
ity 

The luminosity function of galaxies is a tool often 
used to investigate galaxy evolution on cosmological time 
scales. This makes sense as the total luminosity of a galaxy 
seems, also by use of the TF-relation, related to the to- 
tal mass of a galaxy. In determining the local LF one has 
sometimes failed to notice that galaxies are extended and 
that for selection correction one should not treat them 
as point sources. This can lead to a change in the slope 
of the LF as function of redshift and results in the faint 



blue galaxy problem (McGaugh 1994). Another disadvan- 
tage of the one-dimensional LF is that each luminosity 
bin contains galaxies with totally different surface bright- 
nesses and scalelengths (Fig. H). The physical processes in 
disk galaxies seem to be more related to surface bright- 
ness than to total mass (see also Paper IV). Therefore, 
to investigate the distributions related to both the total 
mass and the surface brightness, the luminosity function 
has been divided in several central surface brightness bins 
in Fig. 23. The absolute magnitudes were calculated from 
the values given in Paper I. Figure ^ shows the bivariate 
distribution of a local sample of galaxies and can be useful 
as a reference for high redshift samples observed with the 
Hubble Space Telescope. 



Figure ^ shows that the LF is more or less the same 
for all central surface brightness bins in the B pass- 
band. In the K passband something becomes apparent 
which was already hinted at in the B passband. The LF 
for fainter central surface brightnesses is lower and/or 
shifted to lower absolute luminosities. Figure |2^ is some- 
where in between both options presented in Fig. 3 of Mc- 
Gaugh (1994), which means that both the shape of the 
LF and its normalization depend on the bivariate distri- 
bution of M and /io- Unfortunately, this data is too scarce 
to make a firm quantitative statement, but it is clear that 
further attention should be given to "the LF" . 

At least as important as the number density of the 
local universe is the local luminosity density, which can 
be transformed in a luminous mass density using an M/L 
description. In order to calculate the luminosity density, 
the number density value of each galaxy used in Fig. |2^ 
was given an additional weight depending on its absolute 
luminosity, which results in Fig. |2j. This figure shows the 
total luminosity one expects to find in a random Mpc^ 
from galaxies in the indicated bins of {fiQ^M). The lumi- 
nosities are expressed in solar luminosities per passband, 
calibrated using the absolute solar luminosity values of 
Worthey ( |1994[ ). 

The distribution in the B passband is remarkably flat, 
almost all bins that contain galaxies are equal to within 
1.5 order of magnitude. In the K passband distribution 
there is more structure, as most of the K passband light 
in the local universe comes from higher surface brightness 
galaxies. Due to the scarceness of the data, this figure is 
again more of qualitative than of quantitative interest. 

5.3. Hubble classification 

The Hubble sequence is one of the basic ingredients of 
galaxy formation and evolution schemes, even though the 
underlying physical processes are only partly understood. 
In this section I describe the consequences of some of the 
relations between the structural parameters and Hubble 
type as presented here. 

Figure pO| showed that B/D ratio cannot be used to de- 
termine the Hubble type of face-on systems, which means 
that the classification of edge-on systems is different from 
that of face-on systems. Furthermore, the B/D ratios seem 
to be quite small (< 0.5), even in the K passband where 
the differences in color (and M/L) due to stellar popula- 
tion effects between disk and bulge are minimized. The 
B/D ratios are larger when r^/'' law bulges are used. 
The difference between exponential and r^/'* law bulges 
in terms of generalized exponential profiles is extensively 
discussed in Paper II. Young & Curry ( 1994 ) showed that 
for ellipticals and dwarf ellipticals there is a trend in profile 
shape with luminosity. The brighter ellipticals have r^/* 
like profiles, while fainter (dwarf) ellipticals have more 
exponential like profiles. If this is also true for bulges, 
we might expect early-type spirals to have more centrally 
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Fig. 24. The volume corrected luminosity distribution of galaxies in the (^o,A^)-plane. Each galaxy in Fig. ^ was given an 
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peaked bulge profiles than later type spirals. In Paper II 
it was found that r^^"^ law bulges gave smaller residu- 
als for T<3 than exponential bulges. This picture conflicts 
with the scale independence of Hubble type as seen in the 
K passband data of Fig. |l8|. Within one Hubble type, a 
range in integrated bulge luminosities exists, which should 
result i n diff erent profiles when the model of Young & 
Curry ( 1994 ) is applied to bulges. There is a weak indica- 
tion for this trend, because the galaxies with the brightest 
bulges are slightly better fitted by the r^/^ law bulges 
than by the exponential bulges. The r^/** law bulges never 
give the smallest residuals, not even for the most luminous 
bulges. 



Using an exponential bulge, an important relation is 
found, which is not apparent if an r^/^ law bulge is used: 
h correlates with re, but these parameters do not cor- 
relate with Hubble type (Fig. |l^). This is an important 
relation as it makes the Hubble sequence scale free. Each 
Hubble type comes in a range of different sizes, both in 
terms of diameter and total luminosity. It is an example of 
larger (in the sense of scale size, not mass) galaxies having 
more of everything, larger disk, larger bulge. Due to this 
correlation the scale parameters cancel each other out in 
calculating B/D ratios, which means that a plot of /Xc-/^o 
(oc So/Sqi important in density wave models) versus type 
looks like Fig. ^ with a different scaling. Consequently 
/ie-/io is also a bad diagnostic for Hubble type. 



20 



R.S. de Jong: Near-IR and optical broadband observations of 86 spirals. Ill Statistics of bulge and disk 



The relationship between Hubble type and /io (Fig. |^ 
holds, independent of bulge profile function used. In fact 
the relation holds if one just fits a line to the central re- 
gion of the profile and uses the true central surface bright- 
nesses, because the bulge generally dominates the luminos- 
ity in this central region. One might wonder if in classify- 
ing galaxies one has mainly looked at the surface bright- 
ness of the bulge and not at the B/D ratio. For the earlier 
systems this is harder to accept; the central regions are 
in general overexposed on photographic plates used for 
classifying. It is instructive to know that SO galaxies do 
not fit in this relation. The central surface brightnesses 
of SO galaxies range from ~17.7 to _B-mag arcsec"^, 
estimated from the data of Kormendy (1977) and Peletier 
et al. ( 1990| ). This translates to effective surface bright- 



nesses in the range of 19.5-21 S-mag arcsec" , in accor- 
dance with the two SO galaxies in the current sample, but 
significantly below the trend of the rest of the spiral galax- 
ies (Fig.|lO|). 

If B/D ratio is such a bad diagnostic for Hubble type, 
we are left according to Sandage ( |l96l| ) with only two 
other classification discriminators: 1) the spiral arm struc- 
ture, 2) the pitch angle of the arms. The two remaining 
criteria indicate that the Hubble sequence should be ex- 
plained in terms of spiral arm appearance, even though 
the second criterion might also be in doubt, as measure- 
ments by Kennicutt (1981) showed that pitch angle has 
no tight correlation with morphological type. 

The main theory on spiral structure is the spiral den- 
sity wave theory (Lin & Shu 1964 Roberts et al. 1975 ). 
The fact that Hubble type is a scale free quantity fits into 
this theory. When the bulge and the disk scale with the 
same amount, the shape of the rotation curve stays the 
same, only its amplitude changes. This means that the 
shape of the resonances also scale along with the scale- 
length changes. If bulge and disk scalelengths are corre- 
lated, the shape of the rotation curve is fully determined 
by the relative brightness of the bulge and disk compo- 
nent. Therefore it is harder to understand in the density 
wave theory, why Hubble type does not correlate tightly 
with B/D ratio (Fig. |2^), or to be more precise with (j.^- 
Ho- The strength and the pitch angle of the density wave 
gets modified by the ratio of mass that participates in the 
density wave to the mass that does not. If B/D ratios are 
so small that they hardly could affect the density wave 
(Fig. ^) and if on top of that the B/D ratios and /ig- 
/zq values are only weakly correlated with Hubble type, 
it seems that some modifications to the standard density 
wave model are needed. Maybe a connection between the 
He of the bulge and the distribution of dark matter can 
solve this problem. 



6. Conclusions 

The statistics of the fundamental parameters of 86 spiral 
galaxies have been studied in the optical and the near- 



IR. The use of the near-IR K passband enabled for the 
first time determination of these parameters without being 
hampered by the effects of dust extinction and differences 
in stellar populations. Volume density distributions with 
respect to the fundamental parameters were made, which 
was possible due to the careful selection of the sample. 
The main conclusions are as follows: 

— Freeman's law of a preferred disk central surface 
brightness value needs a modification. Although there 
seems to be a clear upper limit to the central surface 
brightnesses of galaxies, there is no clear limit at the 
faint end of the fiQ distribution. The number of galaxies 
in a volume with a certain /io is only slowly declining 
function of fiQ. 

— The bulge and disk scalelengths are correlated, param- 
eterised by log(rf) = 0.95 log(/i^)- 0.86. This correla- 
tion suggests that the formation of the bulge and the 
disk is coupled. 

— The Hubble classification is related to the surface 
brightness properties of spiral galaxies. However, the 
relations are in general not very tight and can therefore 
not be turned around to give morphological classifica- 
tion. The B/D ratio is not a good indication of Hubble 
type (Fig. po| and ^l|) . The best relation with Hubble 
type found in this study is the one with /ie (Fig. p^ . 
The physical interpretation is difficult, because cause 
and effect are hard to separate. 

— Hubble type is a scale size independent parameter, 
but not a total luminosity independent parameter of 
a galaxy. Therefore, it would be better to divide by 
scale size instead of luminosity to derive scale inde- 
pendent parameters of galaxies in comparisons. The 
suggestion that Hubble type is mainly driven by total 



mass (Zaritsky 1993) seems an oversimplification. To 
carry this point a bit further, it is probably better to 
separate the determination of the LF of galaxies into 
bins which are related to the effective surface bright- 
ness (like Fig. ^3|) than into bins which are related to 
Hubble type. 

Larger samples are needed to enable the parameter- 
ization of the bivariate distributions. There is especially 
need for accurate surface photometry of a large sample of 
galaxies, selected from deep photographic plates providing 
isophotal diameters at very faint levels. 
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